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NOMENCLATURE  -  PART  I 


A 


1  2 
A  («2) 

A.  (ft2) 
a  (ft/sec) 


B 


variable  used  in  variation  of  parameter  solution  to  the  injection 
dynamical  equation 

constant  of  integration 

chamber  cross  section  area,  £  x 

area  of  an  injector  orifice 

local  speed  of  sound 


area  ratio  parameter, 


r  +  l 

2  I2*?  -  *tL 

y + 1/  Ac 


c 

p 


ft-lb 

slug  -  °R 


specific  heat  at  constant  pressure 


Cp£ 

F  (lb) 

G 

gc  (lbm/  slug) 

I  (sec) 
sp 

J 

K 

L  (ft) 

£  (ft) 
c 

£.  (ft) 

K  (ft) 

1 

^^(niole  ) 

M 


injector  orifice  pressure  loss  coefficient 
axial  rocket  engine  thrust 

redefined  in  this  report  as  the  mass  flow  fraction  of  gaseous 
propellant 

unit  conversion  constant  (32.  174) 
theoretical  specific  impulse 

integral  defined  in  the  derivation  of  the  jump  conditions  for  the 
cryogenic  case 

proportionality  constant  in  the  parabolic  approximation  to  P  (r?) 

distance  between  detonation  waves,  277  R/N 

radial  width  of  the  chamber 

characteristic  injector  orifice  length 

width  of  the  nozzle  throat 

molecular  weight 

Mach  number,  v/a 


.  slugs 
m  — — 
sec 


slugs 

sec 


mass  flow  rate 


mass  flow  rate  of  gaseous  propellant 
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NOMENCLATURE  -  PART  I  (continued) 


m 


l£  ' 


slugs1 

sec 


m 


N 

N. 

1 


slugs) 


sec 


mass  flow  rate  of  liquid  propellant 

total  propellant  mass  flow  rate 

number  of  waves  in  the  engine 
number  of  injector  orifices 


P,  Pc  (lbs/ ft  instantaneous  detonation  chamber  pressure 


Pa  (lbs/ft  ) 


atmospheric  pressure 


PaVg  (lbs/ft  )  mean  detonation  chamber  pressure 


PD  (lbs/ft  ) 

Pg  (lbs/ft2) 

AP  (lbs/ft2) 

R  (ft) 

T  (°R) 
t  (sec) 

u,  v,  x 

V  (ft/ sec) 

V  (ft/sec) 
w 

v,  u  (ft/ sec) 

Vi  (ft/ sec) 

V. 

1  avg 
XG 

x,  y  (ft) 

y  (t),  yj(t)  (ft)  circumferential  coordinates  of  the  boundaries  of  a  control 
0  volume  moving  with  the  fluid 


injector  plenum  pressure 
nozzle  exit  pressure 

total  pressure  drop  across  an  injector  orifice 
radius  of  curvature  of  the  chamber  centerline 
absolute  temperature 
time 

dummy  variables  of  integration 

magnitude  of  the  velocity 

absolute  detonation  wave  velocity 

circumferential  and  axial  velocity  components 

injection  velocity  of  the  liquid  propellant 

mean  value  of  V. 

i 

mole  fraction  of  gaseous  propellant 
axial  and  circumferential  coordinates 


a,  13,  w  constants  appearing  in  the  differential  equation  for  0 

y  ratio  of  specific  heats 
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NOMENCLATURE  -  PART  I  (continued) 


?(t)  (ft) 


V 

I 


<P 


Subscripts 

0 

1 

a 

A 

B 

G 

H 

n 

P 

t 

w 

I 


circumferential  coordinates  of  the  detonation  wave 

dimensionless  circumferential  coordinate,  y/L 

dimensionless  axial  coordinate,  x/x 

’  '  n 


mass  density 


nondimensional  perturbation  in  the  liquid  injection  velocity,  V. 


refers  to  condition  at  y  =  0 

refers  to  condition  at  y  =  L 

refers  to  absolute  velocity 

refers  to  the  unburned  propellant 

refers  to  the  burned  propellant 

refers  to  the  gaseous  propellant 

refers  to  the  liquid  propellant 

refers  to  condition  at  beginning  of  nozzle 

refers  to  the  initial  propellant  condition 

refers  to  condition  at  nozzle  throat 

refers  to  the  detonation  wave 

denotes  f/f 

o 


Note:  Absence  of  subscript  refers  to  the  local  condition  in  the  chamber  for  the 
wave  fixed  system  of  coordinates. 
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SUMMARY 


This  report  presents  the  work  accomplished  on  the  rotating  detonation  wave 
engine  feasibility  program  during  the  period  1  March  1963  to  31  May  1963. 

The  work  described  in  this  report  represents  essentially  the  conclusion 
of  the  supporting  experimental  and  theoretical  studies  intended  to  consume 
approximately  the  first  half  of  the  program.  The  last  half  of  the  program  will 
emphasize  the  testing  of  a  nominal  1000  lb-thrust,  rocket  motor  utilizing 
the  detonative  mode  of  combustion  and  employing  hydrogen  and  oxygen  as 
the  propellants. 

The  theoretical  studies  of  an  analytical  model  for  the  internal  gas  dynamics 
associated  with  a  rotating  detonation  wave  engine  have  been  essentially  com¬ 
pleted. 

The  numerical  computations  to  solve  the  differential  equations  for  the 
analytical  model  have  been  carried  out  for  the  no  mixing  case  for  initial 
propellant  temperatures  of  537  R  and  270  R  plots  of  the  various  dimension¬ 
less  gas  dynamic  parameters  versus  the  dimensionless  circumferential 
coordinate,  q,  are  shown.  A  modification  of  the  analytical  model  to 
handle  the  case  of  one  propellant  injected  as  a  liquid  is  presented.  Work 
on  carrying  out  the  necessary  numerical  computations  is  in  progress. 

An  expression  for  the  specific  impulse  of  the  idealized  rotating  detona¬ 
tion  wave  engine  has  been  obtained  Comparison  with  the  ideal  conventional 
rocket  engine  shows  that  the  vacuum  specific  impulse  is  essentially  the 
same  for  both  cases.  A  list  of  sizing  rules  derived  from  the  analytical 
model  are  presented. 

The  experimental  supporting  studies  have  been  essentially  completed. 

Work  is  continuing  on  the  gaseous  100-lb  thrust  motor 

(1)  All  attempts  to  start  the  100-lb  thrust  motor  have  proven  unsuccessful, 
therefore,  a  linear  motor,  similar  in  other  respects  to  the  100-lb 
thrust  motor,  and  having  transparent  combustion  chamber  walls  has 
been  fabricated  to  study  the  problem  systematically. 

(2)  The  experiments  concerning  detonations  in  heterogeneous,  liquid- 
gas  media  have  been  extended  to  include  the  case  of  the  shattering 
of  smaller  (220  -  580  q )  HO  droplets  behind  detonation  waves. 

It  appears  that  these  droplefs  break  up  in  the  same  manner  and 
time  interval  (  ~  10  ju  sec)  as  the  1000  m  droplets  studied  earlier. 

The  appearance  of  a  detached  normal  shock  wave  upstream  of 
the  row  of  shattered  droplets  is  noted  for  this  case. 


xi 


(3)  The  experiments  concerning  the  effect  of  curvature  and  pressure 
relief  has  yielded  two  conclusions.  First,  in  general  it  appears 
that  in  curved  channels  (of  the  same  approximate  dimensions 
as  that  employed  in  the  experiments)  detonation  waves  will  pro¬ 
pagate  at  velocities  (measured  at  the  channel  centerline)  equal 
to  that  observed  in  straight  tubes.  Secondly,  with  no  confine¬ 
ment  (pressure  relief)  on  the  inner  wall,  measured  detonation 
velocities  in  the  curved  channel  appear  to  suffer  a  degradation 
of  no  more  than  7.  5%  compared  to  the  measured  velocities  of 
completely  confined  waves. 


xii 


I.  ANALYTICAL  MODEL  OF  THE  ROTATING 
DETONATION  WAVE  ENGINE 


A.  INTRODUCTION 

In  References  1  and  2  a  simplified  model  of  the  rotating  detonation  wave 
engine  was  developed.  A  numerical  solution  to  the  case  of  complete  instantaneous 
mixing  between  the  burned  and  unburned  propellants  was  given  in  Reference  2. 

In  the  present  report  numerical  solutions  to  the  case  of  no  mixing  between  the 
burned  and  unburned  propellants  are  also  presented. 

In  addition,  the  analytical  model  has  been  modified  to  include  the  case  of  one 
propellant  being  injected  in  the  liquid  state.  The  injection  dynamics  for  a  liquid 
propellant  are  analyzed  and  sizing  rules  for  injector  design  are  derived. 

The  analytical  model  also  has  been  applied  to  the  calculation  of  the  specific 
impulse  of  rotating  detonation  wave  engines.  Comparison  with  the  specific  im¬ 
pulse  of  conventional  engines  indicates  that  an  ideal  rotating  wave  engine  has  the 
same  vacuum  specific  impulse  as  an  equivalent  steady  conventional  rocket  engine, 
but  that  mixing  between  the  burned  and  unburned  propellants  or  inefficiencies  in 
exit  nozzle  performance  would  cause  a  degradation  in  specific  impulse. 

Also,  the  sizing  rules  for  engine  design  derived  from  the  analytical  model 
are  summarized  for  convenient  reference. 

B.  RESULTS  OF  A  NUMERICAL  SOLUTION  TO  THE  NO  MIXING  CASE 

A  digital  computer  program  to  integrate  numerically  the  differential  equa¬ 
tions  fcr  the  no  mixing  case  (see  References  1  and  2)  has  been  developed.  The 
results  fcr  stoichiometric  gaseous  hydrogen- oxygen  injected  at  537°R  and  270°R 
are  presented  m  Figures  1  and  2,  respectively. 

Compared  to  the  complete  mixing  case  the  temperature  of  the  gases  enter¬ 
ing  the  detonation  wave  for  the  no  mixing  case  is  lower  and  represents  the 
ideal  case  of  a  detonation  propagating  in  an  unburned  media.  Note  that  the 
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velocities  v  ,  and  v_  are  measured  relative  to  the  wave  fixed  coordinates. 

A  B 

C.  APPLICATION  OF  THE  ANALYTICAL  MODEL  TO  THE  CRYOGENIC  CASE 

An  important  problem  to  analyze  is  the  situation  where  one  propellant  is 
injected  into  the  detonation  chamber  in  the  liquid  state  (for  example,  liquid  oxygen 
and  gaseous  hydrogen  as  propellants).  For  the  simplified  analysis  given  here, 
only  the  jump  conditions  across  the  detonation  wave  and  the  definition  of  certain 
parameters  need  to  be  modified  in  the  original  model  of  the  no  mixing  case. 

Assumptions 

1.  The  no  mixing  model  is  to  be  used  as  the  basis  of  the  analysis. 

2.  In  control  volume  A  one  propellant  is  in  the  liquid  state,  the  other  is  in  the 
gaseous  state.  The  volume  of  the  liquid  is  neglected  compared  to  the  volume 
of  the  gas. 

3.  There  is  no  interaction  between  the  liquid  drops  and  the  gaseous  propellant, 
i.  e. ,  the  circumferential  coordinate  of  the  drops  remains  stationary  with 
respect  to  the  engine. 

4.  The  liquid  propellant  evaporates  and  reacts  instantaneously  at  the  plane  of 
the  detonation  wave. 

Control  volume  A  will  effectively  contain  only  gaseous  propellant  since  the 
liquid  occupies  negligible  volume  and  is  assumed  not  to  interact  with  the  gaseous 
propellant.  Hence  the  differential  equations  derived  for  the  all-gaseous  case 
will  remain  valid,  except  that  the  parameters  y^,  Cp^,  G,  and  Tp  will  need  to 
be  redefined. 

At  the  detonation  wave  there  is  a  source  of  vaporizing  propellant.  Conse¬ 
quently,  the  hydrodynamic  jump  conditions  across  the  detonation  wave  need 
modification. 
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Detonation  Wave  Jump  Conditions  Including 
a  Source  of  Vaporizing  Propellant 

The  derivation  of  the  jump  conditions  given  here  follows  the  method  presented 
in  Reference  3. 

Consider  a  control  volume  moving  with  the  fluid.  Within  the  control  volume 
there  is  a  detonation  wave  immediately  preceded  by  a  source  of  vaporizing  pro¬ 
pellant.  The  following  is  a  schematic  view  of  this  system  with  respect  to  the 
wall  fixed  coordinates  (see  page  13  of  Reference  1): 
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Propellant  fluid. 


The  conservation  equations  of  mass  and  momentum  for  the  control  volume 
moving  with  the  fluid  are 


Vl) 

d  f  mf 

dt  J  pdy  =  NA; 

yj(t)  c 

yo(l) 

$  J  P  (-\)  “y  =  V  P0 
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Note  that  the  integrals  are  of  the  form 


therefore, 


y0ft) 


J  = 


^  (y,  t)  dy 


yx(t) 


S(t) 


yo(t) 


f  =ft  J  *(y’t)dy  +  l  J  *(y’t)dy 

C(t) 


yj(t) 


y0(t) 


=  J  |^dy  +  £  -  i//(yrt)  yx 


yj(t) 


(+) 


+  ^(y0,  t)  yQ  -  if/  (C  ,  t)  ? 


I  hm\  n  =  ^r,  (V  '  v  )  -  ^  (V  -  V  ) 
(y0  “  yi)  —  °  dt  ^0  w  aQ'  *1  w  ax 


Utilizing  the  velocity  transformation  from  the  wall  fixed  coordinates  to  a 
coordinate  system  moving  with  the  detonation  wave, 

v  =  V  -  v 
w  a 

one  obtains 


lim  dJ 

(yQ  -  yj)  -0  dt  =*0V0~^ 


1V1 
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The  continuity  equation  becomes,  with  =  p, 


pi  vi  +  na"  =  po  vo 


For  the  momentum  equation,  i (/  =  -  p  v  .  Therefore, 

d. 


-  p0  v0  <Vw  -  v0>  +  pl  V1  (Vw  -  vl>  '  P1  -  P0 
Combination  with  the  continuity  equation  yields 


„  2  2 
P,  +p,  V,  +  rr-r—  V  =  P.  +  pn  V- 
111  NA  w  0  ^0  0 

c 


Non-Dimensional  Form  of  the  Jump  Conditions 
In  non-dimensional  form  the  continuity  jump  condition  becomes 


PAX  VAj  +  m£  =  1 


where  m£  is  defined  as 


m* 


PR  v  NA 

B0  Bo  c 


Previously,  the  mass  flow  parameter  G  was  defined  by 


mT 


G  = 


VB0NAc 


where  rhp  is  the  total  propellant  mass  flow  rate.  The  above  relation,  with  G  =  1, 
can  be  viewed  as  a  definition  of  Ac,  since  Afi  in  the  analysis  represents  an  effec¬ 
tive  detonation  wave  area  rather  than  a  chamber  cross-sectional  area.  Also, 


mP  =  mG+m£ 
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Redefine  G  by  the  relation 


m. 


G  = 


PR  VB  NAC 
B0  B0  C 


Hence  one  obtains  the  relationships 


PA  VA  =  G 
A1  A1 


m£  =  !  -  G 


The  non-dimensional  form  of  the  momentum  equation  is 


P1  +  rB  "l  ''a/  *  rB  Vw  =  yB  +  1 


From  the  definition  of  G  this  becomes 


P1  +  yB  G  %  +  >B  (1-G)Vw  =  rB+1 


In  terms  of  mole  fractions  G  is  expressed  as 


G  = 


WGXG 


"^G  XG  +  ^  ^  "  XG) 


where  X„  is  the  mole  fraction  of  gaseous  propellant. 

VJ 

To  solve  for  Pj,  v^,  and  the  integrated  energy  equation  for  control 

volume  A  (see  Reference  1,  page  27)  is  required: 
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Note  that  in  this  case  y A  and  CpA  refer  to  the  specific  heat  ratio  and  specific 
heat  for  the  gaseous  propellant  in  control  volume  A,  G  refers  to  the  mass  fraction 
of  gaseous  propellant,  and  Tp  refers  to  the  injection  temperature  of  the  gaseous 
propellant. 


Combination  of  this  energy  equation  with  the  continuity  and  momentum  jump 
relations  yields  _ _ 


12  2  * 

V*Vr  2  '» 


A  ‘  c  ;j  rA  '  1  y  2 


where  r  is  defined  by 


i  -  rD  [d  -  G)  v  -  l 
B  L  w 


yA+  1 


P,  =  1  -  rB  [G  vAl  +  (1  -  G)  Vw  -  1 1 


yBGyA1 

yAPl 


In  the  numerical  solution  to  this  problem  the  thermochemical  calculations  of 
Reference  4  will  be  used.  Note  that  the  heat  of  vaporization  of  the  liquid  propel¬ 
lant  modifies  the  initial  temperature  used  in  applying  the  thermochemical  results. 

Work  on  a  digital  computer  program  to  solve  these  equations  for  the  cryogenic 
case  is  in  progress. 


D.  INJECTION  DYNAMICS  FOR  A  LIQUID  PROPELLANT 

An  important  aspect  of  engine  design  is  the  dynamics  of  the  injection  of  liquid 
propellants.  It  is  necessary  to  insure  the  injector  design  will  provide  the  desired 
mass  flow  rate  of  liquid  propellant  and  also  distribute  the  liquid  across  the  axial 
dimension  of  the  detonation  chamber.  The  injector  is  assumed  to  consist  of  a 
large  number  of  holes  of  the  form 
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Detonation 

Chamber 


pc(t) 


Following  an  approach  used  by  Morrison,  *  for  steady  flow  through  this  system 
the  pressure  loss  coefficient  Cpf  is  given  by 


AP 


v 


2 


where  AP  is  the  loss  in  total  pressure  from  the  point  where  PD  is  imposed  to  the 
point  where  Pc  exists;  V  is  a  characteristic  velocity.  Choose  the  injection  ve¬ 
locity  V\  to  be  the  characteristic  velocity.  Then  Cpf  is  given  by 


p  _p  _iDv2 
D  C  2  PH  i 

1  W  2 

2  P£  Vi 


for  steady  flow. 

The  momentum  contained  in  the  orifice  at  an  instant  in  time  is  p.  A.  £ .  V.. 

K£  i  i  l 

Hence  the  pressure  head  due  to  the  inertia  of  the  fluid  is  p^  JL  d  V^/dt,  and  the 
pressure  loss  coefficient  for  the  unsteady  flow  case  can  be  written  as 


♦Unpublished  analysis  made  at  this  laboratory  by  R.  B.  Morrison. 
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Pn  -  pr«>  -  *P.  V,2  -  p,  I,  d  V,/dt 


2 

1  Tr  2 

2Pf  Vi 


As  an  approximation  assume  Cp^  is  constant  with  time. 

Note  that  the  only  justification  of  the  assumptions  employed  here  is  that  they 
provide  a  basis  for  obtaining  a  simple  solution  to  the  problem. 


For  the  rotating  detonation  wave  engine,  however,  P^,  fluctuates  very  rapidly 
with  time,  while  may  be  assumed  constant.  (Here  Morrison  assumed  (P^  -  P^,) 
constant. )  The  differential  equation  for  V^t)  is  then 

dV  1  +  cp  P  -  P  (t) 

i  l  ,r  2  D  Cw 

dt  +  2£ .  Vi  p.  £. 

i  £  l 

To  put  this  equation  in  non-dimensional  form  introduce 
PD  =  pr/po’ 


p  =  Vpo’ 


n  =  vwt/L 


Vi  =  Vi/vB  ’ 
i  i  B0 


P£  ~P/PBq’ 


WL 


resulting  in 


dVi+  1  +  \  _  2  ?p-pW 

d”  +  2Vwtj  ‘  ‘vBPjtiVw 


It  has  been  found  that  P(?j)  can  be  closely  approximated  by  a  parabola  of  the 


form 


P(t?)  =  1  -  K  r ,1/2 

where  K  is  a  constant  (e.  g. ,  K  *  1  -  P^). 


V.  (0)  =  V.  (1) 
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The  boundary  condition  is 


Since  the  equation  is  non  linear  an  exact  solution  would  require  numerical 
integration,  for  example  by  the  Runge-Kutta  procedure.  To  obtain  some  insight 
into  the  effects  of  injection  dynamics,  a  perturbation  solution  to  the  differential 
equation  is  obtained. 

Define  0(n)  as  the  fluctuation  in  V.  from  the  mean  value  V. 

'  1  l avg 


<*>  =  V.  -  V. 

1  i  avg 


Assume  0/V.  <<  1.  This  is  equivalent  to  assuming  >>  1. 

The  differential  equation  for  4>  is 


where  w  = 


(1+  Cp.)  V. 
_ i  avg 


d0  ^  0  1/2 

—  +  u xf>  =  a  +  fir] 


c.  v 
vw 


(1  +  Cp.)  V. 
a=s  .  (  1  ayg 


PD^ 


2£iVw  y  B  Pi  h  v 


w 


0  =  — _K_  — 

VB  Pi  h  Vw 


This  differential  equation  can  be  solved  by  the  method  of  variation  of  param¬ 
eters.  The  solution  to  the  homogeneous  equation 


^  +  ux£  =  0 
dr? 


is 


0  =  Ae 


-o crq 


Let  A  =  A(r7)  and  substitute  into  the  non- homogeneous  equation: 
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-o>/7  dA  .  -u .*//  .  -ix’/j 

e  — - a’Ae  +  u-'Ae  -  a  +  fin 

chj  ' 


1/2 


dA  u*/7  _  1/2  u m 

—  =  ae  '  +  fin  e  ' 
d/] 


A  --  —  e^  +  /i  x1,  '  eJ"X  dx  +  A. 

07  j  1 


where  A^  is  a  constant  of  integration.  Hence  4>  is  given  by 


x  a  .  -U7/7  f  1/2  UTX  .  -lx’/7 

0(/j)  =  —  +  fi  e  x  e  dx  +  A^  e 


07 


0 


Aj  is  evaluated  from  the  boundary  condition  fi(0)  =  0(1): 


a  .  a  .  -o7  1/2  uTx  ,  .  -o> 

-  +  A  =  -  +  fi  e  x  e  dx  +  A,  e 

O)  1  u>  1 


0 


A,  =i 


fi 


\ 


1  -i/ 


1/2  UTX  , 
x  e  dx 


Therefore,  the  solution  for  <t>  is 


if  1 

.  .  a  ,  -o)7 1  f  1/2  o/x  ,  /  e  ^  f  1/2  otx 

0(77)  =  -  +  fi  e  l  x  e  dx  +  1 -  x  e 

w  J  /  ,  J 

0  0 


From  the  definition  of  V.  avg, 


dx 
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1 


j  <t>  dr]  =  0 

0 


or, 


-+  p 
w 


I 


-<jjq 


1/2 


wx  , 
e  dx 


dr]  +  — 

'  LO 

we 


1/2  wx 
x  e  dx 


The  second  term  can  be  evaluated  by  applying  the  relation 


v  du 


resulting  in 


a+  |  / 3  =  0 

or, 


(i  +  Cpjv.  2  p^-i  +  Ik 
i  avg  _P _ 3_ 


2  £i  V 


w 


rBPjt  £i  Vw 


V. 

i  avg 


1  +  f  K] 

(i  +  c) 

£ 


Note  that  the  quantity  1  -  2K/3  represents  a  mean  chamber  pressure,  since 
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1 


1 


p 

avg 


j  p 
0 


(1  -Krj1/2)d7  =  l  -f  K 

0 


Hence  K  can  best  be  found  from  K  =  1.  5  (1  -  P  )  where  P  is  the  area  under 
_  avg  avg 

the  P (77)  versus  77  curve  obtained  from  the  engine  analytical  model. 


For  steady  flow  through  the  injector 


p_p  -  —  p  V  2 
D  C  avg  2  i  avg 

1  2 

2  ^i avg 


or, 


2  fPD  "  Pavg] 

'bpi  <1  +  CP(> 

In  the  linearized  case  the  mean  injection  velocity  does  not  differ  from  the 
steady  state  value.  To  judge  the  unsteady  effect  on  V.  a  solution  to  the  non¬ 
linear  differential  equation  should  be  obtained. 


V. 


1  avg 


Sizing  Laws  Derived  from  the  Liquid  Injection  Analysis 
Using  the  derived  relation  for  one  can  express  the  following  sizing 


relations  in  quantitative  form. 
The  total  n 
Therefore, 


1.  The  total  mass  flow  rate  of  liquid  through  the  injector  should  equal  m^. 


m  =  p  N.  A.  V. 
f  £  1  1  1  avg 


where  is  the  number  of  injector  orifices  and  is  the  area  of  each 
orifice.  In  non-dimens io'nal  form  this  becomes 
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or, 


= 


P„  N.  A.  V, 

£  i  t  i  avg 

NA 


^  NiAt  M  lpp  -  Paygi 

mi  -  NAc  V  V1  +  Cp/ 


2.  The  injected  liquid  should  move  a  distance  in  the  same  time  the  deton¬ 
ation  wave  travels  a  distance  L.  Hence 


or, 


V.  V 

1  avg  w 


-  jl  I2  [p°; 

X”  vwT/rBPt(  i  + 


(x  =  x  /L) 
n  n 


E.  THEORETICAL  SPECIFIC  IMPULSE  OF  ROTATING  DETONATION  WAVE 
ENGINES 

To  derive  the  relation  for  the  specific  impulse  of  the  rotating  detonation  wave 
engine  the  well  known  steady  flow  analysis  for  conventional  rocket  engines  is  ap¬ 
plied  to  calculate  the  thrust  of  an  infinitesimal  segment  of  the  engine.  The  result 
is  integrated  over  the  circumference  of  the  engine  and  division  by  the  propellant 
mass  flow  rate  yields  an  expression  for  the  overall  specific  impulse.  It  is  assumed 
that  the  flow  in  the  nozzle  is  expanded  ideally  to  a  given  exit  pressure  P  . 
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Specific  Impulse  of  a  Steady  Flow  Rocket  Engine 

The  specific  impulse,  I  ,  of  a  rocket  engine  is  defined  as  the  thrust  per 

sp 

unit  mass  flow  rate  of  propellant.  It  is  conventional  to  use  pounds-mass  as 
the  unit  of  mass,  then  I  is  given  in  units  of  seconds.  To  make  the  equation 
dimensionally  correct  the  constant  gc  (32. 174  lbm/slug)  is  introduced. 

Consider  a  rocket  engine  mounted  on  a  test  stand: 


Control  Volume  V',  Surface  S 


From  the  momentum  theorem  for  steady  flow  the  rocket  thrust  is  given  by 

F  =£f[pV  (V-  G)  +  Pft]  dS 
S 

where  n  is  the  local  outward  unit  normal  vector  to  the  element  of  surface  area, 
dS.  For  the  control  volume V one  obtains  for  the  magnitude  of  the  thrust 

F  =  p  u  2  A  +(P  -  P  )  A 
e  e  e  e  a  e 

and  since  the  propellant  mass  flow  rate  rhp  equals  pg  ug  Ae>  the  specific  impulse 
is 
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u  (P  -  P  )  A 
t  -.I  /  •  e  e  a  e 

sp  P&c  gc  mpgc 

The  nozzle  exit  velocity  ug  is  found  from  the  energy  equation 

C_,  T  +  J  u  2  =  C„  T 
P  e  2  e  Pc 

and  the  relation  for  isentropic  nozzle  flow, 

y  -  l 


y 

'PC  I 


Hence  the  result  for  u  is 

e 


r 


y  -  l 


2C.T 
P  c 


P  » 

1  -  i  — —  : 

P 

c 


Specific  Impulse  of  a  Rotating  Detonation  Wave  Engine 

To  apply  the  steady  flow  specific  impulse  relation  to  a  rotating  detonation 
wave  engine,  consider  the  flow  field  in  the  wave  fixed  coordinate  system.  The 
flow  is  steady,  then  the  nozzle  flow  can  be  represented  ideally  as  consisting  of 
a  large  number  of  infinitesimal  steady  engines.  Note  that  the  axial  velocity 
component,  u,  has  the  same  value  in  both  the  wall  fixed  and  wave  fixed  coordinate 
systems  and  hence  that  the  thrust  calculation  can  be  performed  in  either  reference 
system.  Consider  one  such  incremental  steady  engine.  The  energy  relation  is 


T 

e 


T 


B 


Note  that  now  T  depends  upon  the  detonation  chamber  circumferential  coordinate 

D 

y.  Also,  for  isentropic  flow  in  the  nozzle, 


T 

e 
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hence 


’'b-1 

P  \  yB 


The  mass  flow  rate  through  an  infinitesimal  nozzle  segment  with  exit  area 


f  dy  is  (see  page  23  of  Reference  1) 


dAt  •  /  2 

N  "  PB  aB  lt  (yB  +  1 


yB  +  1 


Hence  the  thrust  of  this  incremental  nozzle  segment  is 


d'V“VB‘-T| 


rB~  1 
P  \  yB 


PB  aB  *t  y  +  1 


rB*  1 


dy  +  (Pe-Pa)tcdy 


Integrating  over  the  length  L  between  detonation  waves  and  multiplying  by  the 
number  of  waves,  one  obtains  for  the  specific  impulse 


yB  "  1 


=  f 
J 


2  CPB  TB  1  '  [  Pj  PB  aB  ft  >  B  +  1 


VB+  1 

2!^1) 


+  (Pe-Pa)Ody 


Using  the  relations 


m  =PR  v  A  N 

P  B0  B0  C 


-17- 


L  L 


B  = 


yB+1 

'  2<^T) 


yB+1l 


the  non-dimensional  form  of  the  relation  is 


1 

r 


I  =  —  2C„  T„  B  P 

Sp  W  PB  B0  ! 


y  -  1  1 

rB 
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p  p 
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yBSc 
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it  Li 

P  \ 

/P 

6 

Li 
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p  1 

l  c  / 

O  / 

l  a 

There  are  two  special  cases  of  interest: 

1.  Nozzle  exit  pressure  Pg  equals  the  ambient  pressure  P& 


1 

/  1  P  ’ 

I  =  —  /  2  C_  T_  B  P 

sp  gc  V  PB  B0  J 

1  -  [  -1  — 

L  1ppo) 

B 


B 


d  j) 


2.  Ideal  expansion  to  vacuum  conditions  (Pg  =  0) 


I 


sp  vac 


-i  V2CpbTb0BH 


The  analogous  expression  for  a  steady  flow  engine  is 


(I  )  =  —  f 

sp  vac's  gc  Y 


2CP  TC 
c  v  PS  S 
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Hence  the  ratio  of  vacuum  specific  impulse  for  the  rotating  detonation  wave  engine 
to  the  value  for  a  steady  engine  is 


vac^RDWE 
sp  vac^s 


This  ratio  has  been  found  to  be  essentially  unity  for  the  case  of  no  mixing  between 
the  burned  and  unburned  propellant.  However,  there  exist  two  possible  causes  of 
performance  degradation  of  the  rotating  detonation  wave  engine: 

1„  Mixing  between  the  burned  and  unburned  propellant 
2.  Non-ideal  nozzle  performance. 

It  is  anticipated  that  the  extent  of  these  effects  will  be  evaluated  experimentally. 


F.  SUMMARY  OF  ENGINE  SIZING  RULES  DERIVED  FROM  THE  ANALYTICAL 
MODEL—  CRYOGENIC  CASE 

The  following  is  a  summary  of  the  sizing  rules  that  have  been  derived  from 
the  analytical  model  as  applied  to  the  no  mixing,  cryogenic  case: 

1.  Sizing  rules  from  similarity  considerations  (see  Reference  2,  page  16). 


A  *• 
c 


mT 


>BNPo 


rhp  B 


T  2,tRvb  P 


1)Cp  tb 

B  0 


(yB+1 


yB+  1 


Note  that  G  as  defined  in  Reference  2  is  unity  in  this  case.  Also,  there 
are  the  definitions 


L  -  277  R/N 
A  =  x  £ 


c  n  c 
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Note:  Tgp  can  be  found  from  the  results  of  a  thermochemistry  calcula¬ 
tion  for  the  detonation  wave  (see  Reference  4,  for  example)  instead  of 
utilizing  Q  and  Q  as  indicated  in  Reference  2.  This  method  is  also  used 
to  determine  the  molecular  weight  and  specific  heat  ratio  of  the  burned 
propellant. 
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G.  RESULTS  AND  CONCLUSIONS 


The  numerical  calculations  to  solve  the  differential  equations  for  the  rotating 
detonation  wave  engine  analytical  model  as  indicated  in  References  1  and  2  have 
been  carried  out  and  the  results  are  presented  in  Figures  1  and  2  of  Reference  2 
and  in  Figures  1  and  2  of  this  report. 

A  modification  of  the  engine  analytical  model  to  handle  the  case  of  one  propel¬ 
lant  being  injected  as  a  liquid  is  presented  here.  Work  on  carrying  out  the  neces¬ 
sary  numerical  computations  for  this  analysis  is  in  progress. 

An  expression  for  the  specific  impulse  of  the  idealized  rotating  detonation 
wave  engine  has  been  obtained,  and  comparison  with  an  ideal  conventional  rocket 
engine  shows  that  the  theoretical  vacuum  specific  impulse  is  essentially  the  same 
for  both  cases.  More  detailed  comparisons  of  the  relative  performance  of  rotating 
detonation  wave  and  steady  flow  engines  are  being  considered. 

For  convenience,  the  complete  list  of  sizing  rules  derived  from  the  analytical 
model  have  been  summarized. 

Besides  the  remaining  numerical  computations  to  be  completed,  future  work 
concerning  the  analytical  model  chiefly  involves  obtaining  correlations  with  experi¬ 
mental  results. 
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II.  EXPERIMENTAL  STUDIES 


A.  THE  GASEOUS  100-LB  THRUST  MOTOR 

All  attempts  to  start  the  100-lb  thrust  motor  have  proven  unsuccessful. 
Numerous  starting  schemes,  including  the  starting  system  first  used  during  the 
1961  test  series,  have  been  tried.  These  1961  experiments  have  been  duplicated 
as  nearly  as  possible,  even  to  the  point  of  modifying  the  injector  face  to  corres¬ 
pond  with  the  machine  shop  errors  made  on  the  original  injector  face.  Many  of 
the  starting  methods  used  have  resulted  in  a  unidirectional  detonation  wave  that 
quenches  after  the  first  revolution.  Initially,  it  was  believed  that  this  phenomenon 
was  solely  due  to  the  fact  that  the  blockage  in  the  combustion  chamber  was  suffi¬ 
cient  to  stop  the  propagation  of  a  detonation  wave  in  either  direction.  However, 
subsequent  experiments  indicate  that  there  is  heat  addition  occurring  in  the  com¬ 
bustion  chamoer  between  the  first  revolution  of  the  detonation  wave  and  the  closure 
of  the  propellant  valves.  It  is  impossible,  at  this  time,  to  determine  which  of 
these  two  processes  is  the  major  problem. 

In  order  to  solve  the  starting  problem  for  the  100-lb  thrust  motor,  it  is 
imperative  that  the  development  and  propagation  of  a  detonation  wave  be  studied. 
This  study  must  be  made  at  the  same  conditions  that  are  present  in  the  100-lb 
thrust  motor,  and  the  study  can  probaoly  best  be  performed  by  optical  methods. 
Therefore,  the  test  section  shown  in  Figure  3(d)  has  been  fabricated  and  is  now 
ready  for  testing.  The  construction  of  this  linear  motor  is  similar  to  the  100-lb 
thrust  motor,  and  could  be  duplicated,  except  for  the  plexiglas  windows,  by  cut¬ 
ting  the  100-lb  thrust  motor  along  a  radial  line  and  forming  it  into  a  straight  test 
section. 

Installation  of  the  linear  motor  has  been  quite  efficient  because  it  was  designed 
to  utilize  the  same  test  facility  as  the  100-lb  thrust  motor.  The  initial  linear  mo¬ 
tor  test  runs  will  be  in  the  same  manner  as  the  gaseous  100-lb  thrust  motor  runs 
except  the  run  time  will  be  shortened  in  order  to  protect  the  windows. 

The  major  purpose  of  these  tests  is  to  observe  the  development  of  the  detona¬ 
tion  wave  and  its  interaction  with  the  starting  mechanism.  None  of  the  experi¬ 
ments  will  be  concerned  with  the  gas  dynamics  of  the  combustion  chamber  after 
the  detonation  wave  reflects  from  the  test  section  ends,  however,  some  emphasis 
will  be  placed  on  the  injector  face  flow  pattern  prior  to  ignition. 

The  first  starting  arrangement  to  be  tested  will  be  a  spark  plug  igniter  and 
either  a  film  or  brass  shim  stock  blockage  in  the  combustion  chamber.  The 
blockage  will  be  placed  at  various  stations  in  the  comoustion  chamber  with  empha¬ 
sis  placed  on  the  region  where  the  detonation  wave  is  still  in  the  development  stage. 
The  pertinent  data  to  be  observed  during  these  tests  will  be:  interaction  of  the  com¬ 
bustion  zone  with  the  blockage,  duration  of  the  blockage,  and  the  effectiveness  of 
the  blockage.  If  this  scheme  is  successful  in  establishing  a  unidirectional  detona¬ 
tion  wave,  the  blockage  must  also  be  subjected  to  a  fully  developed  detonation  wave 
in  order  to  determine  if  the  starting  mechanism  is  completely  satisfactory. 
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B.  DETONATION  THROUGH  HETEROGENEOUS  LIQUID-GAS  MEDIA 


In  Reference  2  a  series  of  spark- schlieren  photographs  is  presented  showing 
the  breakup  process  of  a  row  of  water  droplets  exposed  to  the  passage  of  a  fully 
developed,  C-J  detonation  wave  under  the  following  conditions: 


Droplet  diameter 
Droplet  spacing 
H2  -  O2  mixture 
Magnification 


1000  m 

3/32  -  1/8  inch 
XHo  *  •  ™ 

2.  £l 


Additional  experimental  information  has  been  obtained  utilizing  the  same 
equipment  and  techniques  concerning  the  breakup  process  involved  with  smaller 
water  droplets  exposed  to  detonation  waves  under  the  following  similar  condi¬ 
tions: 


Mean  droplet  diameter 
Range  in  droplet  diameter 
Droplet  spacing 
h2  -  o2  mixture 
Initial  Pressure 
Initial  Temperature 
Magnification 


400  -  500  p 

220  -  580  p 

3/32  -  1/8  inch 

Xh2  ~  2/3  (stoichiometric) 

One  atmosphere 

Room  temperature 

2.75:1 


The  smaller  droplets  were  produced  by  utilizing  a  smaller  hypodermic  tub¬ 
ing  than  utilized  in  producing  the  1000  p  droplets.  Improvements  in  the  optical 
system  which  included  a  focusing  vernier  ad^'Stment  were  made  in  an  attempt 
to  obtain  more  detail  with  a  greater  magnification..  An  increase  in  magnification 
of  only  10%  was  achieved,  however,  due  primarily  to  two  effects.  It  was  neces¬ 
sary  to  conduct  the  experiments  during  daylight  hours  for  this  study,  consequently, 
excessive  exposure  to  extraneous  light  was  encountered  with  the  extremely  fast 
film  (ASA-3000)  employed.  In  addition,  a  residue  accumulation  on  the  surface 
of  the  test  section  windows  occurred  which  made  it  increasingly  difficult  to  obtain 
sufficient  light  transmission  and  optical  detail.  It  was  essentially  for  these  rea¬ 
sons  that  the  optical  quality  of  the  photographs  was  actually  inferior  to  the  earlier 
photographs  presented  in  Reference  2 


The  series  of  photographs  taken  of  the  row  of  droplets  at  different  time 
intervals  after  the  passage  of  the  wave  yields  the  following  information: 


The  droplets  shatter  significantly  in  the  time  interval  of  from  5  to  10  p- sec 
after  the  passage  of  the  wave  similar  to  that  which  was  observed  for  the  1000  ju 
droplets.  This  conclusion  is  arrived  at  due  to  the  observed  increase  in  the 
apparent  droplet  diameter  by  a  factor  of  three  or  more  in  this  time  interval 
It  is  believed  that  this  apparent  increase  is  due  to  the  shear-type  breakup  process 
observed  by  other  investigators  in  which  the  opaque  zone  is  actually  a  region 
occupied  mainly  by  a  very  dense  population  of  micro-droplets  surrounding  the 
remains  of  the  parent  droplet.  It  is  not  possible  to  establish  conclusively  whether 
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a  portion  of  the  parent  droplet  still  remains  after  the  observed  time  interval  of 
about  10  ii- sec.  What  is  indicated  however,  is  that  the  location  of  the  upstream 
edge  of  the  opaque  region  is  unchanged  from  the  original  position  of  the  row  of 
droplets  up  to  about  10  11- sec  after  the  passage  of  the  wave.  In  the  time  interval 
of  from  10  to  15  ju-sec  after  the  passage  of  the  wave,  however,  a  measurable 
change  in  the  upstream  position  has  occurred  for  the  smaller  droplets  (originally 
~  300  ii  diameter).  It  indicates  that  the  upstream  position  of  the  opaque  zone  has 
accelerated  to  an  average  velocity  in  this  time  interval  of  the  order  of  30  to  60 
ft/ sec.  While  this  is  quite  low  compared  to  the  theoretical  gas  velocity  of  about 
4200  ft/ sec  behind  a  stoichiometric  H£  -  O2  detonation  wave,  it  represents  the 
same  order  of  acceleration  as  that  observed  in  the  shorter  time  interval  (0  to 
10  11- sec)  for  the  downstream  edge  of  the  opaque  zone — a  zone  which  is  believed 
to  be  made  up  of  a  dense  population  of  micro-droplets. 

Figure  4  shows  a  typical  spark  schlieren  photograph  of  the  row  of  droplets 
about  11  ii- sec  after  the  passage  of  the  wave.  The  detonation  wave  is  just  out  of 
the  field  of  view  moving  to  the  right.  For  comparison,  a  photograph  of  a  typical 
row  of  undistrubed  droplets  is  also  shown.  The  closely  spaced  parallel  lines  on 
the  left  side  of  the  photographs  are  reference  wires  spaced  .  078  inches  apart. 

The  weak  normal  shock  wave  just  visible  upstream  of  the  row  of  shattered  drop¬ 
lets  was  not  evident  on  the  earlier  photographs  utilizing  1000  p. droplets.  Several 
possible  reasons  for  the  appearance  of  this  wave  are  under  consideration. 

More  refined  measurements  would  be  required  to  establish  conclusively 
whether  or  not  the  parent  droplets  of  the  order  of  200  to  300  p  initial  diameter 
are  completely  shattered  into  a  spray  of  micro-droplets  in  the  observed  time 
interval.  It  is  believed,  however,  that  more  precise  measurements  of  this  type 
would  yield  significant  information  concerning  the  detailed  droplet  shattering 
processes  occurring  in  an  actual  heterogeneous  detonation  wave,  i.  e. ,  a  detona¬ 
tion  in  which  the  fuel  or  the  oxidizer  component  of  the  reactants  are  originally  in 
the  liquid  droplet  form. 
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C.  GEOMETRICAL  TESTS 

1.  Introduction 

These  experiments  are  essentially  a  conclusion  of  the  work  reported  in 
references  (1)  and  (2)  concerning  the  effect  of  curvature  and  pressure  relief  on 
a  fully  developed  detonation  wave.  The  experiments  utilize  three  geometrical 
configurations: 

(1)  A  curved  channel  with  the  cross-sectional  dimensions 

of  1/2  x  3/8  inch  with  a  radius  of  curvature  of  3.  75  inches 
having  complete  confinement  i.  e.  ,  solid  walls. 

(2)  A  curved  channel  with  the  same  dimensions  as  in  (1) 
above  except  with  the  elimination  of  the  inner  wall 
allowing  a  two-dimensional  expansion  (pressure  relief) 
of  the  burned  gases  behind  the  wave  in  the  inward,  radial 
direction.  Provision  is  made  for  the  placement  of  a  thin 
amyl  acetate  film  along  this  inner  wall  separating  the 
unburned  gases  from  the  atmosphere. 

(3)  The  same  channel  dimensions  as  in  (1)  and  (2)  above 
with  one  window  removed  so  that  the  relief  is  in  the 
axial  direction 

The  majority  of  the  experiments  utilized  a  stoichiometric  mixture  of 
hydrogen  and  oxygen. 

2.  Experimental  Test  Equipment 

A  schematic  diagram  of  the  basic  test  section  with  solid  walls  is  shown 
in  Figure  5  .  A  photograph  of  the  test  section  is  shown  in  Figure  6  .  The 
curved  section  includes  270  .  Provision  is  made  for  the  removal  of  the 
inner  wall  for  the  full  270°. 

A  different  section  including  only  90  0  was  utilized  for  some  of  the  studies 
of  detonations  with  inward,  radial  relief.  This  section  has  provision  for  a 
thin  membrane  (15%  collodion  in  amyl  acetate).  Details  of  the  construction 
and  film  preparation  is  given  in  Reference  2  . 

Each  curved  test  section,  when  used,  is  attached  to  a  straight  detona¬ 
tion  tube  driver  section  of  the  same  cross-sectional  dimensions  as  the 
curved  sections.  The  straight  driver  section  is  nine  feet  long  to  insure  a 
fully  developed  detonation  wave  before  entry  into  the  curved  section.  A 
miniature  glow  plug  is  used  for  ignition  purposes.  A  schematic  diagram 
of  the  basic  system  appears  in  Figure  7  . 
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3.  Experimental  Procedure 

The  experiments  performed  in  the  test  section  with  the  solid  walls  is 
obtained  in  the  following  manner.  A  diaphragm  is  used  to  seal  the  end  of  the 
test  section,  and  the  tube  and  test  section  are  evacuated  by  means  of  a 
vacuum  pump.  Then  the  Hg-  mixture  is  introduced  into  the  tube  and  test 
section  from  a  vessel  containing  the  pre- mixed  combustible  gases.  This 
vessel  is  not  shown  in  Figure  7  .  The  mixture  is  then  ignited. 

The  experiments  performed  in  the  sections  employing  pressure 
relief  are  made  in  either  of  twc  ways 

(1)  The  Hg-  mixture  from  the  vessel  containing 
the  pre-mixed  gases  is  allowed  to  flow  slowly 
through  the  straight  tube  and  test  section  for  at 
least  one  minute  to  purge  the  air  and  other  gases 
completely.  The  valve  to  the  pre-mixed  gas 
vessel  is  closed  and  the  end  of  the  test  section 

is  sealed  immediately  with  masking  tape.  The 
mixture  is  then  ignited. 

(2)  A  flowing  system  of  H~-  O ^  is  introduced  continually 
through  a  mixing  tee  (Figure  7  )  and  metered  by 
means  of  sonic  orifices  in  each  line.  The  mixture 
is  ignited  with  the  system  flowing.  This  method 
was  also  employed  when  utilizing  the  test  section 
with  solid  walls  for  reasons  discussed  later. 

For  each  method,  the  detonation  velocity  is  measured  in  the  straight 
tube  by  means  of  ionization  probes.  In  addition,  the  time  interval  to  the 
spark  light  discharge  is  also  measured.  Two  CMC,  Model  757  BN  time 
interval  counters  are  employed  for  the  two  measurements.  To  photograph 
the  waves,  a  schlieren  system  is  used  employing  six  inch  diameter  mirrors 
instead  of  the  lens  system  shown  (for  simplicity)  in  Figure  7  .  The 
effective  exposure  time  for  the  high-voltage,  capacitor  discharge  light 
source  is  approximately  0,  1  p-sec. 

4.  Results 

Because  of  the  presence  of  small  leaks,  the  final  experiments  uti¬ 
lizing  the  test  section  with  solid  walls  employed  the  flowing  system  for 
consistency  of  results. 

In  addition,  for  the  experiments  exploying  pressure  relief,  the  method 
of  charging  the  tube  by  first  purging  and  then  sealing  the  end  was  found  un¬ 
desirable  because  it  was  determined  that  the  unburned  mixture  partially 
diffuses  through  the  thin  membrane  in  the  time  interval  of  about  three 
seconds  between  the  closing  of  the  valve  and  ignition.  This  effect  was 
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determined  by  noting  that  the  degree  of  wave  curvature  could  be  changed 
by  varying  the  aforementioned  time  interval.  See  Figure  24(a)  of 
Reference  2  which  shows  a  wave  with  relief  showing  this  wave  curva¬ 
ture  effect. 

For  this  reason  the  following  experiments  were  performed  em¬ 
ploying  the  flowing  system.  Figure  8  (a)  shows  a  typical  schlieren 
photograph  of  a  detonation  wave  with  inward  radial  relief  provided 
by  utilizing  a  thin  amyl  acetate  membrane  on  the  inner  wall.  The  wave 
curvature  is  quite  similar  to  a  wave  shown  in  Figure  8(b)  propaga¬ 
ting  in  the  curved  channel  with  solid  walls.  Note  that  the  direction 
of  travel  is  different  for  the  two  waves  i.  e. ,  the  wave  in  Figure  8  (b)  is 
moving  counter-clockwise.  On  the  premise  that  the  degree  of  wave 
curvature  near  the  relief  boundary  is  a  measure  of  the  degree  of  confine¬ 
ment  of  the  wave,  it  is  concluded  that  the  confinement  provided  by  the 
amyl  acetate  membrane  is  essentially  the  same  as  that  provided  by  the 
solid  walls.  It  is  for  this  reason  that  the  experiments  performed  utilizing 
the  membrane  were  not  carried  out  to  conclusion. 

A  composite  distance  versus  relative  time  plot  of  stoichiometric 
H 2-  Og  detonation  waves  in  the  solid  walled  curved  section  employing  a 
flowing  system  is  shown  in  Figure  9  .  The  distance  as  shown  is  obtained 
by  measuring  the  wave  position  along  the  outer  wall.  It  is  apparent  due  to 
the  linear  nature  of  the  average  of  the  experimental  data  that  the  wave 
velocity  is  constant  throughout  the  curved  section.  The  wave  velocity 
obtained  from  the  slope  of  the  average  line  is  9762  ft/sec.  From  similar 
plots  the  velocity  on  the  inner  wall  and  on  the  center  line  are  8541  ft/ sec 
and  9152  ft/sec  respectively.  The  average  wave  velocity  measured  in  the 
straight  tube  before  the  curved  test  section  is  9159  ft/sec,  a  value  very 
close  to  that  observed  along  the  center  line  of  the  curved  section. 

A  schematic  drawing  showing  the  essential  features  of  a  typical 
detonation  wave  in  a  curved  channel  is  shown  in  Figure  10,  obtained 
from  typical  schlieren  photographs  similar  to  Figure  8  (b).  Two  interest¬ 
ing  effects  are  shown  in  addition  to  the  curvature  of  the  wave  from  the 
radial  direction.  First,  it  is  noted  that  stria  exist  in  the  burned  gases 
behind  the  wave  characteristic  of  typical  detonations.  Also  observable 
are  multiply-reflecting  shock  waves  emanating  initially  from  the  outer 
wall  intersection  with  the  detonation  wave  front.  It  is  apparent  that  such 
a  reflected  wave  must  accompany  an  oblique  shock  wave  in  steady  flow. 

Upon  close  scrutinization  of  many  original  photographs  it  appears  that  the" 
observed  stria  indicate  the  approximate  direction  of  the  streamlines  in 
the  burned  gases  because  they  are  deflected  in  the  same  manner  as  that 
predicted  by  two  dimensional  oblique  shock  theory.  It  is  concluded  there¬ 
fore  that  the  observed  shock  system  which  does  not  dissipate  with  the  dis¬ 
tance  travelled  by  the  wave  in  the  curved  tube  is  a  characteristic  of  a  steady 
state  detonation  in  a  curved  channel. 
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It  is  of  interest  to  calculate  the  theoretical  pressure  gradient  in  the  radial 
direction.  This  quantity  is  represented  by  the  following  relation: 


3P2  p2V22 


(1) 


where  pg  and  V-  are  the  pressure,  density  and  absolute  gas  velocity  in  burned 
gases  immediately  behind  a  detonation  wave.  Using  theoretical  values  for  the 
properties  behind  a  stoichiometric,  C-J  detonation  (initially  at  1  atmosphere  and 
room  temperature),  a  value  of  63.  5  psi/ inch  is  obtained  using  relation  (1).  Utiliz¬ 
ing  the  experimental  distance  of  1/2  inch  between  the  outer  and  inner  wall,  a  value 
of  about  32  psi  is  obtained  for  the  overall  radial  pressure  differential  theoretically 
predicted  using  the  1/2  inch  section.  It  is  to  be  noted  that  this  differential  is  about 
12%  of  the  theoretically  predicted  pressure,  Pg,  of  276  psia  immediately  behind  a 
one-dimensional  detonation. 


Another  method  can  also  be  used  tc  predict  the  radial  pressure  gradient  behind 
a  detonation  in  a  curved  channel  by  utilizing  the  momentum  equation  across  the 
wave, 


1 


1 


>1  Mt- 


~r  2  M2 


(2) 


where  and  M„  are  the  Mach  number  in  the  unburned  and  burned  cases  respec¬ 
tively  for  a  one-aimensional  wave.  The  Mach  number  r  ormal  tc  ihe  wave  which 
is  required  by  this  relation,  can  be  obtained  from  measurements  cf  the  wave  angle 
at  points  along  the  wave,  It  was  determined  from  the  experiments  that  the  angu¬ 
lar  deviation  of  the  wave  from  the  radial  direction  varied  from  zero  at  the  inner 
wall  to  16.  4°  at  the  outer  wall.  Using  relation  (2)  above,  the  pressure  differen¬ 
tial  across  the  channel  is  calculated  to  be  about  47  psi  which  compares  favorably 
with  the  value  of  32  psi  predicted  by  Equation  (1). 

Other  tests  were  made  utilizing  the  270°  curved  test  section  with  the  inner 
wall  removed  and  using  a  flowing  stoichiometric  mixture  of  H2-O2.  Figure  11 
is  a  schlieren  photograph  of  a  typical  detonation  wave  in  this  system.  The  wave 
is  approximately  250°  from  the  beginning  of  the  curved  section.  Because  of  the 
lack  of  detail  on  the  photograph,  Figure  12is  included  as  a  sketch  of  the  essential 
details.  With  the  flowing  system  utilized,  a  high  concentration  of  combustible 
mixture  is  maintained  apparently  due  to  centr  ifugal  effects  allowing  the  detonation 
wave  to  be  sustained.  Although  some  mixing  must  be  occurring  between  the  un¬ 
burned  mixture  and  the  atmosphere,  the  effect  must  be  small,  at  least  near  the 
beginning  of  the  curved  section.  Also  iust  visible  in  Figure  llis  the  first  reflected 
wave  originating  from  the  intersection  of  the  detonation  wave  and  the  outer  wall. 
Multiple  reflected  shock  waves  are  not  observed  due  to  the  absence  of  a  reflec¬ 
tion  from  the  gaseous  interface  between  the  burned  gases  and  the  atmosphere. 
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Figure  13shows  a  distance  versus  relative  time  plot  from  a  limited  number 
of  waves.  Again,  in  a  manner  similar  to  the  case  with  the  tube  having  solid  walls, 
the  velocity  (slope)  of  the  average  curve  seems  to  be  linear.  However,  in  compar¬ 
ing  the  detonation  velocity  at  the  channel  centerline  with  the  velocity  obtained  pre¬ 
viously  with  a  solid  wall  inner  boundary,  the  velocity  obtained  with  the  inner  wall 
removed  appears  to  be  about  7.  5%  lower  over  the  complete  length  of  the  curved 
section.  This  is  rather  surprising  inasmuch  as  it  would  appear  that  the  mixing 
effect  between  the  unburned  gases  and  the  atmosphere  would  become  progressively 
more  noticeable  as  the  wave  moves  into  the  curved  section.  From  the  experi¬ 
ments,  apparently  this  effect  is  not  measurable  over  the  distance  involved,  and  the 
degradation  of  about  7.  5%  in  the  detonation  wave  velocity  is  actually  due  primarily 
to  the  relief  effect  of  the  unconfined  gases  behind  the  detonation  wave.  It  is  to  be 
noted  that  certain  theoretical  and  experimental  results  of  a  similar  study*  under¬ 
way  at  this  laboratory  utilizing  detonations  with  relief  in  straight  tubes  tend  to 
support  these  observations. 

Figure  14shows  a  schlieren  photograph  of  a  detonation  wave  propagating  in  a 
channel  with  one  window  removed  so  that  in  effect  an  expansion  of  the  burned  gases 
takes  place  in  the  axial  direction.  The  wave  has  passed  through  the  first  four 
inches  of  the  curved  section.  It  is  apparent  from  the  photograph  that  the  detona¬ 
tion  wave  has  degenerated  into  a  curved  shock  front  followed  by  the  combustion 
zone  (indicated  by  the  turbulent  region  beginning  about  1/4  inch  behind  the  initial 
shock  wave).  It  is  quite  possible  that  the  detonation  wave  is  in  the  process  of 
being  quenched  due  to  the  mixing  processes  which  could  be  more  pronounced  for 
this  case  where  centrifugal  effects  are  not  encountered  in  the  same  stabilizing 
manner  as  in  the  case  with  the  inner  wall  removed. 

5.  Conclusions 

(1)  In  general  it  can  be  stated  that  in  curved  channels  of 
rectangular  crossection  (of  the  same  approximate  dimen¬ 
sions  as  that  utilized  here)  detonation  waves  will  propagate 
at  velocities  (measured  at  the  channel  centerline)  equal 

to  that  observed  in  straight  tubes. 

(2)  With  pressure  relief  on  the  inner  wall,  measured  detona¬ 
tion  velocities  in  the  curved  channel  appear  to  suffer  a 
degradation  of  no  more  than  7.  5%  compared  to  the  meas¬ 
ured  velocities  of  completely  confined  waves.  With  a 
very  thin,  amyl  acetate  membrane  separating  the  unburned 
gases  from  the  atmosphere  along  the  inner  wall  of  a 
curved  channel,  possibly  no  degradation  in  velocity  will  be 
encountered. 


♦Dabora,  E. ,  unpublished  work  in  progress  at  this  laboratory. 
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(3)  Further  detailed  studies  would  certainly  be  warranted  with 
more  carefully  controlled  procedures  and  conditions  than 
utilized  in  these  experiments. 


III.  STUDY  PLANS  FOR  THE  NEXT  QUARTER 


The  theoretical  studies  involving  the  simplified  model  will  be  essentially 
completed  with  the  results  of  the  numerical  computations  of  the  modified 
model  now  in  progress.  The  results  of  the  study  will  be  applied  to  sizing 
considerations  for  the  nominal  1000 -lb  thrust  motor.  The  validity  of  the 
theoretical  model  will  be  examined  when  experimental  results  from 
the  100-lb  thrust  motor  are  available. 

The  work  involving  the  100-lb  thrust  motor  will  be  suspended  until  a 
reasonable  solution  to  the  starting  problem  is  obtained.  The  linear  motor 
with  transparent  walls  will  be  utilized  for  this  purpose.  Following  this, 
a  performance  study  utilizing  the  basic  100-lb  thrust  motor  will  be 
initiated.  Emphasis  will  be  placed  on  the  possibility  of  upgrading  the 
entire  experimental  system  involved  to  obtain  a  nominal  1000-lb  thrust 
from  the  gaseous  motor. 


-31- 


REFERENCES 


1.  Nicholls,  J.  A. ,  and  Cullen,  R.  E. ,  et  al. ,  The  Feasibility  of  a  Rotating 
Detonation  Wave  Rocket  Motor,  Univ.  of  Mich.  Eng.  Res.  Inst. ,  Report 
05179-2-P,  Dec.  1962. 

2.  Nicholls,  J.  A.  ,  and  Cullen,  R.  E. ,  et  al. ,  The  Feasibility  of  a  Rotating 
Detonation  Wave  Rocket  Motor,  Univ.  of  Mich.  Eng.  Res.  Inst.  ,  Report 
05179-3-P,  Mar.  1963. 

3.  Courant,  R. ,  and  Friedricks,  K.  O.  ,  Supersonic  Flow  and  Shock  Waves, 
Inter  science  Publishers,  Inc. ,  New  York,  1948,  pp.  121-124. 

4.  Zeleznik,  F.  J.  ,  and  Gordon,  S.  ,  A  General  IBM  704  or  7090  Computer 
Program  for  Computation  of  Chemical  Equilibrium  Compositions,  Rocket 
Performance,  and  Chapman- Jouguet  Detonations,  NASA  TN  D- 1454,  Oct. 
1962.  (Note:  Some  results  for  stoichiometric  Chapman- Jouguet  detona¬ 
tions  are  presented  in  "Calculation  of  Detonation  Properties  and  Effect 
of  Independent  Parameters  on  Gaseous  Detonations,"  ARS  Journal,  32, 
606  (1962). ) 


Solution  to  the  no  mixing  case’ 


The  dependent  parameters  satis¬ 
fying  both  the  differential  equa¬ 
tions  and  the  hydrodynamic  jump 
conditions  are; 


Cra  ■  0.5497 


Q  =  0.788 
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Figure  1  (a).  Dimensionless  Properties  P,  T^,  Tg,  pB,  and  m  as 
Functions  of  the  Dimensionless  Circumferential  Coordinate,  q,  for 
Tp=  537°R,  Gaseous  No  Mixing  Case. 


Solution  tothe  no  mixing  case 

G  ■  1.0 
TP  •  0.081 
yA  ■  i.4oi 
?B  .  1.15 
Cp.*  0.5497 


The  dependent  parameters  satis¬ 
fying  both  the  differential  equa- 
tionsand  the  hydrodynamic  jump 
conditions  are  * 

B  »  5.76 
Vw  *  1.797 
Q  *  0.788 


Figure  1  (b).  Dimensionless  Properties  M^,  Mg,  v^,  Vg  as  Functions 
of  the  Dimensionless  Circumferential  Coordinate,  r\,  for  Tp=  537°R, 
Gaseous  No  Mixing  Case. 


Solution  to  the  no  mixing  cose 

G  -  1.0 
TP  >  0.04075 
*A  '  15 
iB  •  M5 
Cra»  0.4923 


The  dependent  parameters  satis¬ 
fying  both  the  differential  equa¬ 
tions  and  the  hydrodynamic  jump 
conditions  are- 

B  >  7.608 
VW  •  1,8303 
Q  -0.8037 


Figure  2  (a).  Dimensionless  Properties  P,  TA,  Tg,  pB,  Z,  and  rii  as 
Functions  of  the  Dimensionless  Circumferential  Coordinate,  q, 
for  Tp=  270°R,  Gaseous  No  Mixing  Case. 


Solution  to  th«  no  mixing  case  • 

G  *  1.0 
TP  -0.04075 
h  ■  1.5 
*B  '  "5 
Cpa-  0.4923 


The  dependent  parameters  satis¬ 
fying  both  the  differential  equa- 
tionsand  the  hydrodynamic  jump 
conditions  are  ■ 

B  •  7.608 
Vw  •  1.8303 
Q  •  0.8037 


Figure  2  (b).  Dimensionless  Properties  M.  ,  M  ,  vA,  vn  as  Functions 

A  B  A  B 

of  the  Dimensionless  Circumferential  Coordinate,  r),  for  Tp=  270°R, 

Gaseous  No  Mixing  Case. 
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Figure  3  (a).  Schematic  View  of  the  Linear  Motor. 


Figure  3  (b).  Photograph  of  the  Linear  Motor  on  the  Test  Stand. 


Figure  4.  Schlieren  Photographs  of  the  Shattering  of  H2O  Droplets  Behind 
an  H_  -  0„  Detonation  Wave  (X  =  .  67). 

22  H2 

(a)  Undisturbed  Droplet  Row  (220-580  p  Diameter)  (Retouched) 

(b)  Shattered  Droplets  11-12  p-sec  After  Passage  of  Wave. 


Figure  5.  Schematic  Drawing  of  the  Curved  Detonation  Tube  Section  with 

Solid  Walls. 


Figure  6.  Photograph  of  the  Curved  Detonation  Tube  Section  with  Solid  Walls. 
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Figure  7.  Schematic  Diagram  of  the  Basic  Experimental  System  for 
Studying  Detonation  Waves  in  Curved  Channels. 


Figure  8  (a).  Schlieren  Photograph  of  a  Stoichiometric  H2  -  O2  Detonation 
Wave  in  a  Curved  Channel  Utilizing  a  Thin,  Amvl-Acetate  Membrane  for 
Inward  Radial  Relief  Employing  a  Flowing  System. 


Figure  8  (b).  Schlieren  Photograph  of  a  Stoichiometric  H2  -  O2  Detonation 
Wave  in  a  Curved  Channel  with  Complete  Solid  Wall  Confinement  Employ¬ 
ing  a  Flowing  System. 
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Figure  10.  Interpretive  Sketch  of  a  Stoichiometric  H2  -  O2  Detonation  Wave 
in  a  Curved  Channel  with  Complete  Solid  Wall  Confinement. 


Figure  11.  Schlieren  Photograph  of  a  Stoichiometric  H2  -  O2  Detonation 
Wave  in  a  Curved  Channel  with  the  Inner  Wall  Removed  Employing  a  Flow¬ 
ing  System. 
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Figure  12.  Interpretive  Sketch  of  a  Stoichiometric  H2  -  O2  Detonation  Wave 
in  a  Curved  Channel  with  the  Inner  Wall  Removed  Employing  a  Flowing  System. 
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Figure  14.  Schlieren  Photograph  of  a  Stoichiometric  H9  -  O2  Detonation 
Wave  in  a  Curved  Channel  with  One  Window  Removed  (Relief  in  the  Axial 

Direction) 


